The Solar Cycle: A new prediction technique based on 

logarithmic values 



O 

(N 
> 

Q 

o 

(N 

6 



Z. L. Du 

Key Laboratory of Solar Activity, National Astronomical Observatories, Chinese Academy of Sciences, 

Beijing 100012, China 
zlduSnao . cas . en 



ABSTRACT 

A new prediction technique based on logarithmic values is proposed to predict the maximum 
amplitude {R-m) of a solar cycle from the preceding minimum aa geomagnetic index (aoniin)- 
The correlation between In/Jm and Inaomin [r — 0.92) is slightly stronger than that between 
i?m and aomin = 0.90). From this method, cycle 24 is predicted to have a peak size of 
i?m(24) = 81.7(1 ± 13.2%). If the suggested error in aa (3 nT) before 1957 is corrected, the 
correlation coefficient between R^a and aa^i^ [r = 0.94) will be slightly higher, and the peak of 
cycle 24 is predicted much lower, i?in(24) = 52.5 ± 13.1. Therefore, the prediction of i?„i based 
on the relationship between R^a and aamin depends greatly on the accurate measurement of aa. 
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1. Introduction 

Predicting the strength of an upcoming solar 
cycle (i?m) is important in both solar physics and 
space weather. A variety of methods have been 
used to do so, of which some are based on statis- 
tics and some others are related to phys ics (iKane 
2007t ICameron fc Schussleil l2007t IPesn cll 2008 



Hiremat hI 12008: iTlatovl l2009t iMesserotti et al 



2009 : iPetrovavl [2010t IDu fc Wand I2OI0I ). A re- 
liable prediction of R^ m ay test model s for ex- 
plaining the solar cycle ( Pesnell 2008f ). Var- 
ious s olar dynamo models ( e.g 
20061 : 



Dikpati et al 



Choudhuri et al. 20071) have been pro- 



posed to explain the solar cycle but the pre- 
dictiv e skill of R^ needs to be checked infu- 
tur e (ICameron fc Schiis"slMl 120071 : IPesnelll l2008t 
D^l2011ah . Based on the Sola r Dyna mo Ampli- 
tude (SODA) index, ISchattenI (|2005l ) predicted 
that the peak sunspot numb er of the current cycle 



(24) will be low, at - 80. 'Dikp ati et all (|2006l ) 



predicted that the peak size of cycle 24 will be 
30 %- 50% higher than that of cycle 23 based on 
a mo dified flux-transpo r t dyn amo model. In con- 
trast, IChoudhuri et al.l (|2007l ) predicted that the 



peak size of cycle 24 will be 30 %- 50% lower than 
that of cycle 23 based on a flux-transport dynamo 
model. 

bmce [Ohl (Il966[l found a high correlation 
between the minimum aa geomagnetic activity 
(aamin) in the declining phase of a solar cycle and 
the maximum sunspot number of the succeeding 
cycle (i?m), a great many papers related to this 



(Brown & Williams 


1969; Kane' 


20101 Wilson 


,1990: .Hathaway & Wilson,2006;.Charvatoval2009l: 


IWane & SheelevI 200! 


The level of geomagnetic 



activity near the time of solar activity minimum 
has been shown to be a good indicator for the 
ampli tude of the following solar activity maxi- 



mum (I0hllll976': 'Wilson''l990: Laydcn et al. T991| 



Thompsonlll993.: ,Hathawav fc Wilson .2006; .Kan 



20101 ). This method based on a solar dynamo 
concept that the geomagnetic activity during the 
declining phase of the preceding cycle or at the 
cycle minimum provides approximately a measure 
of the poloidal solar magnetic field that generates 



the to roidal field for the next cycle (jSchatten et al 
19781) . 

When geomagnetic precursor methods are ap- 



1 



plied to the current cycle (24), some discrepancies 
are sh own for different authors. iHathawav &: Wilson 
(|2006l ) predicted i?m(24) = 160 ± 25 using the 
I component of aa by subtract ing the linear R 
component with from aa (jFeynman Il982h . 
Dabas et al. I (l2008l) employed the number of geo- 
magnetic disturbed days prior to the minimum 
of the sunsp ot cycle, and p r edict -Rm (24) 



Wang fc SheelevI (jiooi) predicted 



124 ± 23. 

i?m(24) — 97 ± 25 based on the total open flux 
at sunspot minimum, which is derived from the 
historical aa index by removing the contribution 
of the solar wind speed. 

The correlation coefficients between and 
the geomagnetic-based paramete r s are usually 
very high from 0.8 dohll Il976t IWilsonI Il990l: 
Thompson 1993; 'Shastri' 'l998'; 'Kane '2010") up to 
0.97 [Laydcn et al. .1991: LaMOAfc Richa rd 199j 
Hathawav et al.lTl999l : lOabas et alTl2008l ). How- 
ever, a high correlation coefficient does not al- 
ways yield an ac curate predi ction, such as in the 
case of cycle 23 (|Kanel 120071 ). It was found that 
the correlation coefficient between and aomin 
varies roughly in a cycle of about 44-year and 
that the prediction error based on this method 
when the correlation coefficient decreases is much 
larger tha n that when the correlati on coefficient 
increases (|Du et al.ll2009t lDull2011al) . 

Conventionally, the correlation between i?ni 
and aOmin is analyzed by a linear relationship 
(Section [5]) , and cycle 19 is viewed as anomalous 
or an 'outlier' due to the very great R^ ever seen. 
However, by analyzing the relationship between 
the logarithms of i?,n and aa^m in Section |3l cycle 
19 is no longer anomalous from the scatter points 
of lni?i„ versus InaOmin- Whether correcting the 
suggested error (3 nT) in aa before 1957 has great 
influences on the prediction of i?„i based on the 
relationship between and aa^nin (Section [1]). 
The results are briefly discussed and summarized 
in Section [5] 

2. Linear relationship between R^ and 

This study uses the annual values of geomag- 
netic aa index computed from the 3-hourly K 
indices at two near-antipodal midlatitude sta- 



tions (|Mavaudlll973 iLovdlioilh since 186^1 and 
the equivalent ones from measurements taken in 
Finla n d from 1844 to 18 67 (jNevanlinna fc Kataia 
I993I : iNevanlinna 2004), and the annual values 
of the International sunspot number (i?z) since 
184^1 produced by the Solar Influences Data Anal- 
ysis Center (SIDC), World Data Center for the 
Sunspot Index, at the Royal Observatory of Bel- 
gium. The maximum amplitude of sunspot cycle 
(i?ni) and the preceding aa minimum (aa^in) are 
listed in Table [1] 

Conveniently, one employs the linear relation- 
ship between R^ and aamin to predict the former. 
Figure [T] depicts the scatter plot between i?m and 
aftmin (triangles). The dotted line indicates the 
linear fit of Rm to aamm, 



Rm = 12.9 ± 14.7 + (7.84 ± 1.06)aan 



(1) 



where the values following ± indicate the stan- 
dard deviation. The correlation coefficient be- 
tween i?ni and aflinin is very high, r = 0.90 at the 
99% level of confidence. From this equation and 
aamin(24) = 8.7, the peak size of cycle 24 is pre- 
dicted to be i?m(24) = 81.2 ±16.2, where a = 16.2 
is the standard deviation of fitting, defined by 



N-1 



(2) 



where R[ is the fitted value of i?m by Equation ([T]) 
and TV = 15 is the number of data pairs. 

It is seen in Fig. [l]that the point of cycle 19 is 
far above the fitting line, Ai?ni(19) = 42.0. There- 
fore, cycle 19 is often called as an 'outlier' ( Kand 
'2OO7') and this point may, occasionally, be deleted 
in order to obtain a better correlation. 

3. Relationship between lni?m and Inaamin 

Now, we analyze the scatter plot between the 
logarithms of R^ and aamin, as shown in Fig. [5] 
(triangles) . 

The correlation coefficient of lni?m with Inaamin 
is r = 0.92 at the 99% level of confidence, slightly 
stronger than that of i?ni with aamin (0.90). The 



iftp://ftp.ngdc.noaa.gov/STP/SOLAR.DATA/RELATED- 
JNDICES/AAJNDEX/ 
^http: / /www. sidc.be/sunspot-data/ 
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Fig. 1. — Scatter plot of i?in vs. aomin (triangles) 
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Fig. 2. — Scatter plot of Ini?,,! vs. Inaomin (trian- 
gles). 

least-squares-fit regression equation (dotted line) 
is 

lni?f = 2.58 ± 0.26 + (0.84 ± 0.10)lnoOmi„, (3) 

where lni?f denotes the fitted value of lni?ni. The 
standard deviation of fitting is erf = 0.132. This 
equation is equivalent to the form of a power-law, 



p „ „2.58„„0.84 



(4) 



implying that i?ni does not depend completely lin- 
early on aflniin. The error of lni?t is 



e = Alni?f = lni?f - \nR„ 



(5) 



The values of Rf and e are listed in Table [T] 

One can see from Table [1] that the maximum 
relative error occurs in cycle 19, |e(19)| = 26.8%, 
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"Prom Fig.[2]and Eqs. Q and JSj. 
'From Fig. [3] and Eqs. lO and Q. 



which is only slightly larger than that in cycle 15, 
|£(15)| = 25.0%. Therefore, cycle 19 seems to be 
not an 'outlie r' as cycle 15 in view of the relative 
error. In fact, Ramesh fc Lakshmil ( 2011 ) proved, 
through a thorough analysis of the linear relation- 
ship between R^ and the preceding sunspot min- 
imum (iimin), that cycle 19 is not an outlier — it 
is more appropriate to be called as an anomalous. 

From Equation ([3]), the peak sunspot number 
for cycle 24 can be predicted: IniJf (24) = 4.403 ± 
0.132 or i?f(24) = 81.7(1 ± 13.2%), close to that 
by the linear relationship (81.2). 



4. Using the corrected aa 

The aa index was suggested to exist an er- 



1957 ( 


Nevanlinna & Kataia 


1993: Lukianova et al. 


20091; 


Svakaard et al.l 12004 


). In this section, we 



corrected the suggested error in aa by adding 3 nT 
to aOmin for cyclcs 9-19 and re-examine the pre- 
vious results. The Scatter plot between i?m and 
the corrected aoniin is shown in Fig. [3] (triangles). 

The correlation coefficient between and 
aajjiin is now r — 0.94 at the 99% level of con- 
fidence, slightly stronger than that using the un- 
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Fig. 3. — Scatter plot of i?,n vs. the corrected 
oamin (triangles) by adding 3 nT to aa before 1957 



corrected aamin (0.90). The linear regression equa- 
tion (dotted line) is 

i?,„ -30.5 ± 15.8 + (9.54 ± 1.00)aamin, (6) 

and the standard deviation of fitting is cr = 13.1. 
Table [T] lists the fitted result (i?p) and the predic- 
tion error of Rp, 



ARp — Rp — Rn 



(7) 



From Equation ([6]), the peak sunspot number 
for the next cycle (24) is predicted to be i?,n(24) = 
52.5 ± 13.1, much lower than that using the un- 
corrected aamin in Section [2] (81.2 ± 16.2). There- 
fore, the prediction of R^n based on the relation- 
ship between R^ and aa^nin depends greatly on 
the accurate measurement of aa, that is, whether 
correcting the suggested error in aa before 1957. 

5. Discussions and Conclusions 

Studying the variations in the 11-yr solar cy- 
cle may help to understand the f ormation and 
dynamo mech anism of the cycle (Parker 1955t 
Babcocklll961l ). It has long been noted that the 
11-yr Schwabe cycle is close to the synodic pe- 
riod of t he co-alignments of the Earth, Ve nus, and 
Jupiter (|Woodl Il975t [Crandpierrd Il996l ). How- 
ever, it is uncertainty whether the planetary 
tidal force can tri gger the dynamo mechanism 
( Grandpierrd 1996[) as the acceleration due to 
planetary tidal force is much smaller than the 



observed acceleration at the level of tachocline 
( de Jaeer fc Versteeghl [20o3). In the dynamo 
mechanism, the differential rotation in the so- 
lar convective envelope transforms the poloidal 
magnetic field structure into toroidal magnetic 
field structure which leads to the format i on of 
sunspots due to Coriolis force ( Parkeil 19551 



Babcockl ll96l': 'Schatten et al.ll l978l : iDikpati et al 
20061 : IChoudhuri et al. 2007). Dynamo mod 



els can reproduce certain features of the cycle 
(e.g., sunspot butterfly diagrams), but the pre- 
dictive skill of Rm has not been checked so far 



(ICameron fc Schiissleil 120071 : IPesnelll 120081: |d 
2011al ). As the actual observational time series 



of poloidal field (a vailable only since the mid- 
1970s) is not long (jChoudhuri et al.l l2007l ). the 
geomagnetic activity around the cycle minimum 
is used as a mea sure to estimate the p ol oidal solar 
magnetic field ( Schatten et al. 19781 ). Ijavaraiahl 



(|2008[ ) found that R^ is well correlated with the 
sum of the sunspot group areas in the 0° - 10° 
latitude interval both of the Sun's northern hemi- 
sphere near the minimum of the previous cycle 
(r = 0.95) and of the southern hemisphere just af- 
ter the time of the maximum of the previous cycle 
cycle (r = 0.97). Recently, iTlatovl (|2009f ) sug- 
gested that the parameter G = S(l/A^g)^, defined 
by the number of sunspot groups Ng > 1, may 
be useful for calibration of the residual magnetic 
poloidal fields, as the amplitude of G is highly 
correlated with i?m at one and a half solar cycles 
later (r = 0.96). 

Conventionally, the relationship between Rm 
and aamin is analyzed linearly. The upcoming 
Rm is predicted by extrapolating the linear regres- 
sion equation from a least-squares-fit algorithm 
(i?f), and its uncertainty is estimated by the stan- 
dard deviation (ct) as the actual prediction error 
(Ai?f = R{ — i?m) has not been known until the 
cycle is over. Thus, the prediction is usually ex- 
pressed in the form of i?ni = Rpicr (2cr), regarding 
the 68% (95%) level of confidence. 

Giving the uncertainty in an absolute measure 
(cr) is enough in most circumstances. If a predic- 
tion error (Ai?ni) is less than 20, this prediction is 
usually thought as a successful one. In some cases, 
however, it may alternately be better to describe 
the uncertainty in a relative form. Suppose that 
the prediction errors are the same for two predic- 
tions (i?mi, i?m2), Ai?m = 20. If i?mi > 100, 
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this prediction is rather successful as its relative 
prediction error is less than 20%. However, if 
-Rm2 < 50, this prediction is terrible as its relative 
prediction error is larger than 40%. Therefore, 
showing the prediction error in a relative form is a 
better choice, particularly in comparison with two 
or more predictions. 

This study examined the relationship between 
lni?m and Inaamin, with a correlation coeffi- 
cient (r — 0.92) slightly higher than that for 
the linear relationship between and aamin 
(r = 0.90). The standard deviation of fit- 
ting so obtained (at) refers directly to the rela- 
tive standard deviation. From this method, the 
peak sunspot number for cycle 24 is predicted 
to be i?m(24) = 81.7(1 ± 13.2%), near to that 
from a modified Gaussian function (72±11, 'Du| 
201 Idl). that from the sun spot minimum (85±17, 



Ramesh fc Lakshmil 2011 ) , and that from the sum 
of the sunspot group areas in the 0° - 10° latitude 



interv al of the previous cycle (87 ± 7, iJavaraiah 
20081 ). In fact, t he logarithmic was often used 



the studies of Waldmeierl (1933). 



If the suggested error in gg (3 nT) before 1957 



(iNevanlinna & Kata 


all993; Lukianova et al.ll2009l 


Svakaard et al. 2004 


) is corrected, the correlation 



coefficient between R^^ and ggmin (j — 0.94) will 
be slightly higher than that using the un-corrected 
aamin {r — 0.90). From this method, the peak 
sunspot number for the next cycle (24) is pre- 
dicted to be i?m(24) = 52.5 ± 13.1. It is close 
to that f rom long-t erm trends of sunspot activ- 
ity (55.5, Tan l201l[ ). lower than both that by an 



autoregressive model (110 ± 11, iHirematll 2008 ) 



and that by the G parameter (135 ± 12. iTlatov 



2009t) . The prediction (52.5) is lower than that us 



ing the to tal open fliix derived fr om the gg index 
(97 ± 25. 1 Wang fc SheelevI [iooil ) . that using the 



number of g,eomag netic disturbed days (124 ± 23 
IPabas et al.l l2008[ i. a nd that using the I compo 



nent of aa (160 ± 25, iHathawav fc WilsonI [20061) . 
It should be pointed out that the prediction using 
the co rrected aa ( 52.5). which is close to that by 
Kand (58.0 ± 25.0. I2010I ). is much lower than that 
using the uncorrected gg (81.2 ± 16.2). There- 
fore, the accurate measurement of gg is crucial to 
predict Rm when using the relationship between 
i?ni and the preceding ggmin- Whether correcting 
the suggested error in aa before 1957 may lead 
to great discrepancies in the prediction of R^ by 



using the above relationship. 

Accurately predicting the peak size of a upcom- 
ing sunspot cycle is a difficulty task as the monthly 
sunspot numbers may have systematic u ncertain- 



ties of about 25% (jVitinskii et al.lll986l) . In fact 



the relationship between aa and i?z is very com- 
plex and the current aa yalue ma y be related to the 
past solar activities (lDull2011bl l3). refiect ing long- 
term evolu tion characteristics of the Sun's m ag- 
netic field (jLockwood et al.lll999l ; lTlatovll2009[) . 
Main conclusions can be drawn as follows. 

1. The correlation between lni?m and Inaamin 
(r = 0.92) is slightly stronger than that be- 
tween i?m and aainin('' — 0.90). From this 
method, the i?ni for cycle 24 is predicted to 
be i?m(24) = 81.7(1 ± 13.2%). 



2. The prediction of i?m based on the relation- 
ship between R^ and ggmin depends greatly 
on the accurate measurement of gg. If the 
suggested error in aa (3 nT) before 1957 is 
corrected, the correlation coefficient between 
i?ni and aamin {r = 0.94) will be slightly 
higher, and the peak size of cycle 24 will be 
predicted much lower, i?m(24) = 52. 5± 13.1. 
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